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Abstract | Few-layer black phosphorus was recently rediscovered as a narrow-bandgap atomically 
thin semiconductor and has already attracted unprecedented attention due to its interesting proper¬ 
ties. One feature of this material that sets it apart from other atomically thin crystals is its structural 
in-plane anisotropy which manifests in strongly anisotropic transport characteristics. However, tra¬ 
ditional angle-resolved conductance measurements present a challenge for nanoscale systems such 
as black phosphorus, calling for new approaches in precision studies of transport anisotropy. Here 
we show that the nonlocal response, being exponentially sensitive to the anisotropy value, provides 
a powerful tool for determining the anisotropy. This is established by combining measurements of 
the orient at ion-dependent nonlocal resistance response with the analysis based on the anamorphosis 
relations. We demonstrate that the nonlocal response can differ by orders of magnitude for different 
crystallographic directions even when the anisotropy is at most order-one, allowing us to extract 
accurate anisotropy values. 
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Black phosphorus (BP), a layered crystal allotrope of 
phosphorus, is composed of corrugated sheets of phospho¬ 
rene held together by weak van der Waals interactions.^ 
In phosphorene sheets, each phosphorus atom covalently 
bonds to its three nearest neighbors, forming warped 
hexagons in a non-planar, chair conformation. This sp^- 
type bonding results in a distinctly anisotropic crystal 
structure: phosphorene adopts an accordion-like folded 
hexagonal lattice with grooves running along the y axis 
(with zigzag termination). Such crystal structure leads 
to highly anisotropic energy bandselectrons and holes 
of phosphorene are nearly an order of magnitude heav¬ 
ier in the y direction (zigzag) as compared to the x 
direction (armchair). This results in an interesting 
in-plane anisotropy in electronic, optical and acoustic 
properties.^ Further, electrical conductivity and carrier 
mobility are higher along the x axis than along the y 
axis, whereas the opposite anisotropy sign is predicted 
for thermal conductance, a property that can potentially 
lead to high-efficiency thermoelectric devices.^ The in¬ 
plane anisotropy is also of interest for plasmonics, since 
the plasmon resonance frequency could be tuned with the 
light polarization due to the highly anisotropic nature of 
plasmons in BP.^ However, while being of high interest 
because of their promise for device research and engi¬ 
neering, the anisotropic properties are currently lacking 
reliable characterization techniques. Here we develop a 
new approach that is particularly well suited for nanosys¬ 
tems where other known techniques can perform poorly 
or lack precision. 

Transport anisotropy was studied recently in a two- 
terminal configuration by angle-resolved conductance 
measurements.^’^ Refs. 6 and 7 employed a circular geom¬ 
etry, measuring two-terminal conductance between sets 
of opposite contacts spaced by 30° or 45°. The angular 
dependence for the two-terminal conductance was fitted 


to an expression (j{0) = (Jxx cos^ 0 + ayy sin^ 0, with 0 
the angle between the chosen pair of opposite contacts 
and the x axis. From the measured dependence, both 
the crystallographic orientation and the values axx^ o-yy 
could be inferred, giving the anisotropy value 

^ ^ XX ! ^yy ( 1 ) 

ranging from 1.5 to 1.8. Unfortunately, such direct two- 
terminal measurements are prone to inaccuracy due to 
contact resistance and variability in device geometry, and 
thus are expected to yield relatively large errors in the 
inferred anisotropy values. 

Here we employ a different approach. We obtain the in¬ 
plane anisotropy from nonlocal measurements in a four- 
terminal geometry, in which a pair of current leads is 
positioned far from the pair of voltage probes. Our ra¬ 
tionale for resorting to such a method is the exponential 
sensitivity of the nonlocal signal to the anisotropy value. 
The exponential character of the nonlocal response fol¬ 
lows from the seminal van der Pauw theorem,^ which 
stipulates a relation between four-terminal resistance val¬ 
ues Rij^ki — Vij/Iki obtained for the contacts k, I and ij 
used as current leads and voltage probes, respectively. 
This theorem connects the nonlocal resistance values for 
different combinations of voltage and current contacts 
in arbitrary-shape homogeneous 2D samples. For an 
isotropic conductor of the sheet conductivity a it reads 
g-7rcri?i2,34 g-7rcrit:23,4i — where the values Rij^ki de¬ 
pend on contact placement. Applied to long samples of 
dimensions L x W, L/W ^ 1 with two pairs of contacts 
placed at opposite ends, the van der Pauw relations pre¬ 
dict an exponential dependence of the nonlocal response 
on the aspect ratio 

Rni oc , p = a~^. (2) 

Below we extend this relation to anisotropic media using 
the anamorphosis (reshaping) theorem (see Appendix). 
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The anamorphosis theorem stipulates that the multi¬ 
terminal resistance values Rij^ki remain unchanged upon 
reshaping the sample via an anisotropic rescaling of the 
coordinates along the principal conductivity axes^”^^ 

X = 7 ^/^a;', y = 7 = , (3) 

and simultaneously letting the conductivity of the re¬ 
shaped system to be isotropic and equal 

O-' = {(7xX<:ryyY/‘^. (4) 

For long samples with the major axis parallel to the x 
or y axis the reshaping gives L/W = /W' ^ 1. 

The van der Pauw relation then predicts an exponential 
nonlocal response altered by the anisotropy as 


ocp'exp(^-:^^ , (xp'exp(^-^^ (5) 

where we defined p' = 1/cr'. Here and R^^^ denote 
the nonlocal response for long samples aligned with the 
X OT y axis, respectively. Our argument based on reshap¬ 
ing indicates robustness of the exponential 7 dependence. 
The order-one prefactors in Eq.(5), which in general are 
not universal, will be analyzed below for a strip geometry. 

The strong dependence on 7 suggests that nonlocal 
measurements afford a precision tool for determining the 
anisotropy. To exploit this opportunity we fabricated 
multiterminal Hall bars aligned with the x and y crys¬ 
tallographic axes of black phosphorous, see Fig. la-c. 
As discussed below, conductance anisotropy can then be 
extracted from the ratio of the nonlocal and local resis¬ 
tance. We found this approach to be reliable and robust 
despite a large variability in the nonlocal response val¬ 
ues that could differ by orders of magnitude for different 
crystallographic directions. 

To better understand the impact of anisotropy on the 
nonlocal response, Eq.(5), we first consider an isotropic 
conductor shaped into a Hall bar as in Eig. lb. When cur¬ 
rent is passed between, say, contacts 2 and 3, it mostly 
flows straight across the device, with an exponentially 
small fraction flowing sideways, as the continuity equa¬ 
tion requires that current does not drop to zero abruptly 
in the device channel. This produces stray currents which 
give rise to a voltage drop across remote contacts 4 and 5 
(Eig. lb). Analysis predicts a nonlocal resistance which 
decays exponentially with the separation L between cur¬ 
rent leads and voltage probes: 


p — ^4,5 2/9 

rinl = — - = -In COth 

h,3 TT 




where the asymptotic exponential form gives accuracy 
better than 0.1% already for L = W. 

A more complex behavior arises for an anisotropic con¬ 
ductor, where the anisotropic conductivity defines an 
easy (high-conductivity) and a hard (low-conductivity) 
direction for current flow. In this case, current applied 



FIG. 1. Fabrication and characterization of crystal-axes- 
aligned devices of hBN-encapsulated few-layer phosphorene. 
a,b) Device micrograph and schematics: two orthogonal Hall 
bars oriented along the x (easy) and y (hard) crystal axes, c) 
Device cross-section: BP, encapsulated with hBN on SiO^/Si 
substrate, connected to gold contacts through tunnel barriers 
provided by the monolayer hBN. d) Polarization-dependent 
Raman spectra measured along two principal axes of hBN- 
encapsulated BP (x axis, red; y axis, black). Inset: Angle- 
resolved Raman signature / B 2 g used to determine the crys¬ 
tallographic orientation of a BP flake, e) Four-terminal con¬ 
ductivity of hBN-encapsulated BP vs. gate voltage measured 
at T = 5K and bias voltage of O.IV {x axis, red; y axis, blue). 
Inset: temperature dependence of the field-effect mobility for 
electrons (solid squares) and holes (empty squares). 


along an easy direction is “pinched” into a relatively nar¬ 
row path, which suppresses stray currents and reduces 
nonlocal response. In contrast, the current applied along 
a hard direction tends to spread out due to an “easy” 
transverse flow short-circuiting the direct “hard” flow. In 
this case, enhanced stray currents lead to a large nonlocal 
response. The resulting suppression and enhancement of 
i?ni in the two regimes is exponential^ Eq.(5), with the 
corresponding exponents differing by the anisotropy fac¬ 
tor A = 7 ^ (see Appendix). 

Next, we describe fabrication of the BP devices used 
in this work. Since few-layer black phosphorus is known 
to lack stability under ambient conditions,we exfoli¬ 
ated BP crystals (purchased from HQ Graphene) in the 
inert atmosphere of a glove box. In order to protect BP 
from oxygen and moisture, crystals were encapsulated 
with hexagonal boron nitride (hBN), using the “dry peel” 
technique detailed elsewhere.Starting with freshly 
exfoliated BP on a polymer film, a thin (« 10 nm thick¬ 
ness) crystal was identified and picked up with a mono- 
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layer hBN on a polymer membrane.The resulting stack 
was then transferred onto a 30-80 nm hBN which served 
as an atomically flat substrate. The dry transfer method 
ensures an extremely clean interface between BP and 
hBN crystals and protects the former from degradation 
when air-sensitive few-layer BP is exposed to ambient 
conditions for many months^^. To avoid the degrada¬ 
tion of BP surface during the fabrication procedures, the 
micromechanical exfoliation and transfer processes were 
performed in an argon environment provided by a glove 
box with levels of H 2 O and O 2 below 0.1 ppm. The final 
assembly on 290 nm SiOx/p-Si wafer was taken out of the 
glove box for subsequent microfabrication processing. 

We used polarization-dependent Raman spectroscopy 
to establish the BP crystal axes orientation. Measure¬ 
ments were performed at room temperature in air using 
a confocal Witec spectrometer. Laser spot size of 0.5 /im 
was obtained by using a Nikon lOOx objective. Samples 
were illuminated with green (2.41 eV) laser at 100 /iW 
intensity. A high-density grating (1800 g/mm) was used 
to achieve energy resolution better than 2 cm“^. Angu¬ 
lar control of the sample, relative to the laser polariza¬ 
tion, was made with a rotating stage giving an uncer¬ 
tainty of about ±2 degrees. Black phosphorus features 
three main Raman modes Ag for the 2 ; axis, out-of¬ 
plane phonons, for in-plane vibrations along the x 
axis (armchair), and B 2 g for in-plane vibrations along 
the y axis (zigzag). The amplitude of the two in-plane 
modes is polarization-dependent with intensity maxi¬ 
mal for laser polarization along the x axis (armchair), see 
Fig. Id. From the angular dependence of the normalized 
mode intensity A^ / B 2 g we extracted the crystallographic 
orientation of BP crystals, as shown in Fig. Id inset. 

Subsequently, multiterminal Hall bar devices were de¬ 
fined via standard microfabrication techniques - e-beam 
lithography and metal deposition followed by reactive ion 
etching - with the device channel aligned with the x or 
y axis. Fig. la-c. Because the encapsulating hBN mono- 
layer has relatively low (< 1 kO-/im^) tunnel resistivity^^, 
the electrodes (5nm Cr adhesion layer followed by 100 
nm An) were deposited directly onto it. Fig. Ic. Two 
devices (A and B) were studied in this work, each device 
comprising two orthogonal Hall bars defined along the x 
and y crystallographic axes. Devices A and B had chan¬ 
nel length-to-width ratio L/W of 3 and 2.5, respectively, 
with W = 1 /im for both devices. Here we define the 
length L as the distance between voltage probes (e.g. 2 
and 4 or 8 and 10 in Fig.lb), the same length scale will be 
relevant for nonlocal measurements; see below. To char¬ 
acterize transport properties we measured four-terminal 
conductivity along both the x and y axes. DC conduc¬ 
tivity was obtained by passing current through contact 
1 to contact 11 and measuring voltage drop across, for 
example, contacts 2 and 4 for cjxx and contacts 8 and 10 
for (Tyy^ see Fig. lb. Conductivity values were calculated 

as (Jxx — ^XX L/W , (Jyy = Ryy L/W . 

Our field-effect transistor devices showed ambipolar 
behavior with somewhat enhanced hole-type transport 










FIG. 2. Transport anisotropy of black phosphorus devices, 
T = 5K. a,b) Resistance vs. gate voltage measured for 
samples A and B along principal crystallographic axes (red: 
X axis, blue: y axis, see panel c inset). Solid lines: lo¬ 
cal resistivity pxx and pyy. Dashed lines: nonlocal resis¬ 
tance Rn\. Noise level 0.05D. c) The apparent anisotropy 
Aapparent = cTxx/cryy measured as the ratio of the local (longi¬ 
tudinal) conductivities of the Hall bars aligned with the x and 
y axes, d) The anisotropy coefficient (averaged for Vg vary¬ 
ing from —20V to —GOV) determined from the ratio of the 
nonlocal to local resistance, using piocai = Pxx for and 
Piocai = Pyy for R///'• luset shows the anisotropy coefficient 
vs. gate voltage for both samples and both crystallographic 
axes (sample A: x axis, red, y axis, blue; sample B: x axis, 
magenta, y axis, black). 


(a likely result of p-type doping during crystal growth), 
as illustrated in Fig. le for device B. The overall de¬ 
vice performance was comparable with that reported 
previously. We extracted field-effect mobilities pxxi Ryy 
from the slopes 5(7/5Vg in Fig. le as /i = C~^5g j5Vg. 
Here C is the device channel-to-gate capacitance C = 
eeo/d^ where d is the dielectric thickness. The capaci¬ 
tance was determined independently using Hall measure¬ 
ments. The field-effect mobilities plotted vs. tempera¬ 
ture (see Fig. le inset) are nearly constant below 80K, 
decreasing as T~^ above lOOK with a varying from 0.6 
to 1.5 depending on the carrier type and crystallographic 
direction. We interpret the T~^ as the onset of electron- 
phonon scattering. 

At low temperatures (5-50K) the hole mobility val¬ 
ues were approximately 1500 and 500cm^V“^s“^ for 
transport along the x and y axes, respectively. The 
corresponding electron mobility values were somewhat 
smaller, approximately 1200 and 200 cm^V“^s“^, respec¬ 
tively. The anisotropy of mobility was therefore between 
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2 and 4. Accounting for the Drude relation between 
mobility and conductivity, a = ne/r, yields the ratio 
A = (7^^jGyy as 1.5 ± 0.1 and 2.6 ± 0.1 for devices A and 
B, respectively (see Fig. 2c). Such a large spread in the 
apparent anisotropy values illustrates the pitfalls of this 
direct approach for accurate determination of anisotropy, 
Eq.(l). The large spread arises even when the different 
devices used are made from the same crystal, and with a 
well-defined Hall bar geometry. 

One can identify two main reasons for the direct 
method to lack precision^^. First, the value A is ob¬ 
tained from conductivity measurements on physically dif¬ 
ferent samples with somewhat different conductivity ten¬ 
sors. Second, even for well-defined geometry the aspect 
ratio might vary slightly from sample to sample. On the 
contrary, with our new method, both local and nonlo¬ 
cal measurements are done on the same device, thus we 
eliminate these two sources of inaccuracy. Further, the 
exponential sensitivity of the nonlocal response to the 
anisotropy value gives the nonlocal method a distinct 
advantage over local measurements in finding the true 
conductance anisotropy. 

Next, we proceed to discuss the nonlocal resistance 
measurements. Due to the high contact resistance which 
is at least an order of magnitude larger than the channel 
resistance, we are able to reliably assess nonlocal resis¬ 
tances only in a region of Vg from -lOV to -GOV, where 
the contact resistance is relatively small. We obtain non¬ 
local resistance by passing the current transverse the de¬ 
vice channel (e.g., from contact 2 to 3 in Fig. lb) and by 
measuring a voltage drop across the nearest pair of poten¬ 
tial contacts (e.g., between contacts 4 and 5 in Fig. lb). 
Figs. 2a and 2b summarize results of our measurements 
for both local and nonlocal resistances using the two 
studied devices and both principal crystallographic axes. 
Nonlocal resistance was several orders of magnitude lower 
than the local resistance, and was also considerably lower 
along the y axis than along the x axis. Comparing panels 
a and b in Fig. 2, the difference in the nonlocal resistance 
values measured along the x and y axes is larger for the 
device with LjW = ‘^ than for the one with L/W = 2.5. 

This behavior can be exemplified by four-terminal 
measurements in device B (in our discussion below we 
refer to the schematics shown in Fig. la,b). Nonlocal 
resistance measured along the x axis, R^i = ^ 4 , 5 / 12 , 3 , 
is a factor of 5-20 larger than the one measured along 
the y axis, i?ni = V 9 po/^ 7,8 5 see dashed red and blue 
lines in Fig. 2b. Indeed, as discussed above, since cur¬ 
rent between contacts 2 and 3 flows along the hard (low- 
conductivity) direction it tends to spread along the Hall 
bar channel (x-axis), resulting in a larger voltage drop 
across contacts 4 and 5. In contrast, current between 
contacts 7 and 8 flows along the easy (high-conductivity) 
direction, and is therefore pinched into a relatively nar¬ 
row channel and spreads much less in the hard direction 
along the Hall bar channel (y-axis). This results in a 
much smaller voltage drop at remote contacts 9 and 10. 

As discussed above, an anisotropic conductor can be 


treated as an isotropic one after an appropriate coor¬ 
dinate transformation (see Appendix). This is accom¬ 
plished by stretching/shrinking the device dimensions by 
a factor of 7 = (Jxxj^^yy and replacing the anisotropic 
conductivity tensor by an “average conductivity value”, 
Eq.(4), giving 
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( 7 ) 

respectively. Hence, for a fixed geometry, the ratio of 
the nonlocal and local resistance depends only on the 
anisotropy coefficient 7 . This enables us to extract 7 
directly by combining the local and nonlocal measure¬ 
ments on a single Hall bar. Notably, Eig. 2d inset shows 
that 7 values are very similar for both devices and both 
crystallographic orientations of the Hall bars, and almost 
independent of gate voltage. Small deviations in 7 for 
Vg > —20V stem from complications of current pathways 
as the Eermi level approaches band gap. Eig. 2d summa¬ 
rizes data from Eigs. 2a and b and demonstrates the close 
agreement of anisotropy coefficient for both studied de¬ 
vices yielding 7 = (Jxxj^yy ~ 1-29 ± 0.05 and, hence, 
the anisotropy value A = o^xl^yy ~ 1.66 ± 0.1. 

It is interesting to compare these results with the 
anisotropy in bulk BP. Anisotropy of Hall mobility mea¬ 
sured at about 50K for p-type crystals was about 3.3 
(Ref. 22 ), which is notably larger than the values found 
in our experiment as well as the values reported in Refs. 6 
and 7. The twice smaller anisotropy of in-plane conduc¬ 
tance of a few-layer phosphorene as compared to the bulk 
BP can be attributed to the surface scattering, which 
reduces the anisotropy of electron-phonon interaction.^^ 
The anisotropic electron-phonon interactions as well as 
the anisotropic properties of acoustic phonons in black 
phosphorus could also lead to the anisotropy in Drude 
relaxation times Tyglr^x^^ The latter could explain the 
difference between the anisotropies of conductance and 
effective mass in BP. Erom previous measurements^^ and 
theory^^’^^ the ratio of effective masses along the x and y 
directions was estimated at 7.7 ± 1.2, which is 

much larger than the conductance anisotropy measured 
in this work. At constant carrier concentrations, the dif¬ 
ference of the two anisotropies could be accounted for 
by the anisotropy Tyg/T^x ~ 4.4 ± 1 . 1 . It would be in¬ 
teresting to find out how these anisotropies evolve with 
decreasing the number of layers down to monolayer phos¬ 
phorene but, unfortunately, such encapsulated BP de¬ 
vices will require remote-controlled microfabrication in 
an inert atmosphere^^ which, at present, appears to be 
very challenging. 

Lastly, we note that the nonlocal resistance observed 
in BP has purely classical explanation, thus differing it 
from nonlocal responses observed in systems with non¬ 
trivial band properties due to either spin-orbit coupling 
(spin Hall effect) or non-zero Berry curvature (anoma¬ 
lous/valley Hall effect). In the latter two cases, the non¬ 
local resistance is a manifestation of a combination of di- 
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rect and inverse spin/anomalous Hall effects, induced by 
the presence of spin- or valley-polarized currents. These 
effects should play only a minor role (if any) in the nonlo¬ 
cal transport in BP because phosphorus is a light element 
and the spin-orbit terms do not contribute appreciably 
to BP band structure.Moreover, topology of a band 
structure near the P point of BP Brillouin zone is triv¬ 
ial (single valley),leading to zero Berry curvature and, 
thus, to the absence of topological currents. 

Summing up, we have described a new approach to 
characterize transport anisotropy of layered nanomateri¬ 
als which has distinct advantages over previously devel¬ 
oped approaches. Our method leverages on novel mea¬ 
surement schemes (nonlocal response) through which the 
effect of anisotropy can be sharply enhanced and accessed 
directly without resorting to traditional angle-resolved 
approaches. Namely, our method utilizes the combina¬ 
tion of local and nonlocal transport measurements and 
relies on the exponential sensitivity of nonlocal resistance 
to the anisotropy coefficient. Using this method we find 
that the anisotropy of in-plane conductance Oxxl^yy for 
thin BP crystals amounts to ~ 1.66 ± 0.1. This method 
can easily be adapted to other potentially anisotropic 2D 
crystals (e.g. ReSe 2 and ReS 2 ) as well as for 2D materi¬ 
als and heterostructures where in-plane anisotropy might 
be induced via strain or charge density waves. 
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Appendix A: The Anamorphosis Theorem 

Here we derive a simple and general relation that con¬ 
nects multi-terminal resistance values for an anisotropic 
conductor and those for a reshaped isotropic conductor of 
conductivity a' = {c^xx^^yyY^"^• This relation, thereafter 
referred to as the anamorphosis (reshaping) theorem^ is 
valid for homogeneous samples of an arbitrary shape with 
the contacts placed in an arbitrary fashion, either at the 
sample boundary or in the interior. The theorem stipu¬ 
lates that the multiterminal resistance values Rij^ku de¬ 
fined as the potential difference on the contacts i and j 
induced by a unit current fed through contacts k and /, 
remains unchanged upon sample reshaping: 

^ij,kl = Rij,kl‘ (^ 1 ) 


Here Rij^ki = Vij/hr R'ij^ki = ^Ij/^ki where the un¬ 
primed and primed quantities correspond to the original 
and reshaped sample, respectively. Eq.(Al) provides a 
simple and model-independent relation between different 
transport measurement results. Unlike the van der Pauw 
theorem^ which connects multiterminal resistance values 
obtained for the same sample by a permutation of cur¬ 
rent and voltage contacts, the result in Eq.(Al) relates 
the values obtained for different samples albeit with no 
contact permutation. 

To establish the result in Eq.(Al), we recall that the 
conductivity of an anisotropic 2D conductor, in the ab¬ 
sence of magnetic field, is described by a second rank 
tensor. Without loss of generality, we can pick the axes 
of our coordinate system to be aligned with the principal 
conductivity axes, writing 


In these coordinates, the equations describing ohmic 
transport read 

jx ~ ^xxRxi jy — ^yyRy 

where E = (Ex^Ey) = —V(j). Here the current j, elec¬ 
tric field E, and the potential (j) spatial dependence are 
obtained from a solution of a boundary value problem 
posed through the current continuity relation 

diji = 0. (A4) 

This equation, which holds inside the sample material, 
must be solved together with the boundary condition 
specifying fixed potential values at the contacts 0^, i = 
1...4 and the tangential-current condition 1 x j = 0 at the 
nonconducting parts of the boundary (here 1 is a vector 
tangential to the boundary). 

In order to map the anisotropic problem to an isotropic 
one we introduce an anisotropic scaling transformation 

x = 7]x', y = n~^y', r] = (A5) 

We note that under this transformation the continuity 
relation given in Eq.(A4) continues to hold provided the 
current components are redefined as 

3x = jy = V~^fy (A6) 

Crucially, j' satisfies the current-field relations for an 
isotropic conductor of conductivity value a': 

J-; = a'E',, j'y = cj'E'y, E' = (E;, E'y) = -V'0 (A7) 

where V' = {ri~^dx'^gdy'). 

Next, it is straightforward to check that the reshaping 
transformation in Eq.(A5) leaves the tangential-current 
boundary condition 1 x j = 0 unchanged. Indeed, upon 
rescaling in Eq.(A5) the tangential vector changes as 
= {h^x^h~^^y)^ yielding 1' x j' = 0 (our tangential vec¬ 
tors 1 and 1' are not normalized.) We therefore conclude 
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that new transport equations are identical to those for 
an isotropic material of conductivity a'. 

Lastly, we note that the definitions of voltages and cur¬ 
rents at the contacts remain unchanged upon rescaling: 
VI = Vk and 

I'^ = J dV xy = J d\xj = Ik, (A 8 ) 
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where the integral is taken along the boundary of the kth 
contact before and after rescaling. As a result, the multi¬ 
terminal resistance values for the transformed geometry, 
in which both the sample and the contacts are reshaped, 
are identical to those in the original anisotropic sample: 


V ' V- • 

TDf _ _ D 

Tf T ^ 


ij,kl5 


which proves the reshaping theorem, Eq.(Al). 

As a sanity check, consider a conducting strip of width 
W parallel to the x axis which carries constant uni¬ 
form current /. Voltage drop between probes on the 
strip edge separated by a distance L is then given by 
V = {L/Waxx)I- After the reshaping transformation, 
Eq.(A5), the new strip width and the new distance be¬ 
tween voltage probes are W' = r]W and L' = r]~^L. Ac¬ 
counting for the isotropic conductivity a' of the reshaped 
conductor we find that voltage remains unchanged upon 
reshaping: 


V' = = 

W'a' 


Wax 


-/ = V, 


(A9) 


which is in full accord with our anamorphosis theorem. 

Next, we consider currents and potentials for a nonlo¬ 
cal measurement in the same strip. The corresponding 
contact placement is illustrated in Eig.3. Performing the 
reshaping transformation to a strip with an isotropic con¬ 
ductivity a', as above, we arrive at a problem which has 
been solved elsewhere.In this case the nonlocal voltage 
induced at a distance L' from current leads is given by 



TT 


/ COshTT^ + 1 
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where p' = 1/a'. This relation can also be written as 

(All) 


Eor L' W' this gives 

(A 12 ) 

TT 

These expressions, according to the relation (Al) 
proven above, also describe the nonlocal response in 
an anisotropic system, provided that the corresponding 
lengthscales are related via W' = pW, L' = p~^L. This 
provides a derivation of the expression for nonlocal resis¬ 
tance used in the main text. 


FIG. 3. Schematics of anisotropic conductor modelled as in¬ 
finite strip of width W. Principal axes of conductivity tensor 
(cTxcc, ayy) are aligned with the x and y axes. Current I is ap¬ 
plied between source (S) and drain (D) electrodes. Voltage 
probes Vi and V 2 are placed at a distance L away from S — D 
pair. 


Appendix B: Direct derivation of the nonlocal 
response for a strip geometry 

Here we derive the nonlocal response for a strip geom¬ 
etry by a direct method that does not rely on the reshap¬ 
ing theorem. We proceed by combining the current-field 
relations, Eq.(A3), with the continuity equation diji = 0 , 
which gives an anisotropic Laplace’s equation for the po¬ 
tential 


^^ j ° 

This equation must be solved in a strip 

- 00 < X < 00 , -yo <y <yo, yo = W/2 (B2) 

(see Eig.3) with the source and drain current leads po¬ 
sitioned at (x^y) = ( 0 , ±^o) and potential probes lo¬ 
cated at a distance L on the opposite edges of the strip, 
{x,y) = {L,±yo). 

Accounting for translational invariance of the strip 
(B2), we perform Eourier transform via ^{k^y) = 

/ dxe~'^^^(j){x^y)^ converting Eq.(Bl) into a a second- 
order ordinary differential equation 

^k, y) = 0, K = 7 fc, (B3) 

where 7 = {axx /as above. Eq.(B3) has a general 
solution 


^{k,y) = Aue-y ^ Bue-^y. (B4) 


Current entering and leaving the system through the 
source and drain leads is modeled by a Dirac delta func¬ 
tion in the boundary conditions for the current density 
normal component, jy{x) = 16{x) at y = ±yo. Writing 
these relations as d(f)k/dy = —I/ayy and solving for 
Bk, gives 


^k = -Bk 


I 


(B5) 


2K,ayy cosh hiyo 
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Plugging these values in Eq.(B4) we find 


is identical to that in Eq.(AlO). Eor L W this gives 




I sinh hiy 
ka' cosh Ki/o 


(B 6 ) 


Rn\ 



(Bll) 


where we used the relation i^cryy = ka'. Inverting Eourier 
transform yields a solution to the original problem 


7^ dk e sinh hzy 


(B7) 


Eor the potential difference across the strip at a distance 
X = L from current leads, V{L) = (f{L, —yo) — ^{L^yo), 
this gives 


V{L) 



k 


(B 8 ) 


This expression can be evaluated with the help of a con¬ 
tour integral 




(B9) 

where a > 0 and the integration path runs along — oo < 
z < + 00 , returning back at infinity. Applied to V{L), it 
yields 


y(L) = 5 (BIO) 

which is in agreement with the result noted above, 
Eq.(All). The above expression for nonlocal response 


Erom this result we obtain the ratio of nonlocal resistance 
i?nb measured for a strip aligned with the x {y) axis, and 
the corresponding resistivity pxx {pyy)'- 


R 


(x) 

nl 


Pxx 


47 f ttL \ 


R. 


(y) 

nl 


= — exp 

Pyy \ W 

(B12) 

With the quantities on the left hand side known from 
measurements, each of these nonlinear equations can be 
solved individually to obtain the 7 value. This can be 
conveniently written in terms of the so-called Lambert 
Wn function: 


7rI/7 



/ 4Z/ Pxx \ 



pyy \ 


(B13) 

The function Wo is defined as the inverse of the func¬ 
tion f{x) = xe^ in the domain x > 0 (the subscript 
0 denotes the principal real branch of the function). 
The Lambert function is easily accessible from most of 
mathematics software products, including OriginPro^^, 
MATLAB^^ and Wolfram Mathematical^. This proce¬ 
dure gives two independent estimates of the anisotropy 
value A = axxl^yy Notably, despite the fact that in our 

measurements and R^^^ typically differ by orders of 
magnitude, the relations (B13) yield nearly identical 7 
values for the Hall bars oriented along different crystal 
axes on the same device. Similar values are also found 
for measurements on different devices. 
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